One-pion production in neutrino-nucleus collisions 
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O ' Abstract. We use our model for neutrino pion production on the nucleon to study pion production on a nucleus. The model 

^rfH ^ is conveniently modified to include in-medium corrections and its validity is extended up to 2GeV neutrino energies by 

^^ , the inclusion of new resonant contributions in the production process. Our results are compared with recent MiniBooNE data 

^^ ' measured in mineral oil. Our total cross sections are below data for neutrino energies above ~ 1 GeV. As with other theoretical 

calculations, the agreement with data improves if we neglect pion final state interaction. This is also the case for differential 

cross sections convoluted over the neutrino flux. 
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The MiniBooNE Collaboration has recently published one pion production cross sections on mineral oil {CH2) 

^ ' by Vfi/v^ neutrinos with energies below 2 GeV [1, 2, 3]. These are the first pion production cross sections to 

be measured since the old bubble chamber experiments carried out at Argonne National Laboratory (ANL) [4, 5] 

and Brookhaven National Laboratory (BNL) [6]. The latter were measured in deuterium where nuclear effects are 

^^ I small [7, 8]. MiniBooNE data poses an extra problem to theoretical models due to the expected relevance of in- 

medium modifications and final state interaction (FSI) in carbon. In Ref. [9], the use of a formation time/zone, that 
reduces the impact of FSI, leads to a good agreement with the shape of different neutral current (NC) tt" production 

^^ , differential cross sections. Charged current (CC) single pion production off '^C for neutrino energies up to 1 GeV is 

analyzed in Ref. [10]. Their results for total cross sections are below MiniBooNE data in the high neutrino energy 
region (0.8 — IGeV) and the agreement improves if FSI is neglected. A different approach valid only in the low 
Q^ region is presented in Ref. [11]. There the authors evaluate pion production by neutrinos in the low Q^ region 

^^ I and for neutrinos in the energy range 0.5 ~ 2 GeV. The model is based on partial conservation of the axial current 

H ■ (PCAC) hypothesis, the conserved vector current (CVC) hypothesis and the use of experimental cross section data at 

. . .' the nucleon level. The agreement found with MiniBooNE data is good for Q^ values up to 0.2 GeV^. In Ref. [12] the 

authors use the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model finding that total cross sections measured by 
MiniBooNE are higher than theoretical ones for neutrino energies above 0.8 ^ 0.9 GeV. As in Ref. [10], the agreement 
with data for total and different flux averaged differential cross sections is better if pion FSI is neglected. However, 
as also shown in Ref. [12], the same FSI model applied to pion photoproduction is able to give a fair reproduction of 
experiment in that case. 

In this contribution we address the problem of pion production in a nucleus starting from our pion production model 
at the nucleon level taken from Refs. [13, 8]. In order to better compare to MiniBooNE data we extend the model up 
to 2 GeV neutrino energies, well above the A resonance region for which it was originally developed. Above the Delta 
region also the Z)i3(1520) resonance plays a role [14] and in the present calculation we include its contribution. We 
also take into account in-medium corrections to the production process. Those are Pauli-blocking and Fermi motion 
and the important corrections that originate from A resonance modification inside the nuclear medium. Another issue 
is pion FSI for which we use a simulation program that follows the work done in Ref. [15] where a general simulation 
code for inclusive pion nucleus reactions was developed. In some of the channels coherent pion production is also 
possible and to evaluate its contribution we shall take our results in Ref. [16] that uses the model we derived in 
Ref. [17]. 



PION PRODUCTION MODEL 



Our full model for one pion production on the nucleon is depicted in Fig. 1 . It contains the dominant A-pole resonance 
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FIGURE 1. Model for the W^N -^ N'n reaction. We have Direct and crossed A(1232)— and nucleon pole terms, contact and 
pion pole contribution, and the pion-in-flight term [13]. In this case we also include direct and crossed Di3-pole terms. 

term (direct and crossed) and background terms required by chiral symmetry. On top of that we add now the direct 
and crossed D^-pole terms. The background terms are the leading contributions of a SU{2) nonlinear a model 
supplemented with well known form factors in a way that respects both conservation of vector current and the partial 
conservation of axial current hypotheses. All the details on the A and background terms can be found in Refs. [13, 8]. 
In Ref. [8] we followed the work in Ref. [18] and we made a combined fit of the dominant nucleon-to-Delta axial 
form factor to ANL and BNL data including both full deuteron effects and flux normalization uncertainties. In Fig. 2 
we show the results of that fit compared to experimental data. The axial nucleon-to-Delta form factors obtained in 
Ref. [8] are the ones we use in the present calculation. As for the D13 resonance contribution, all details will be given 
elsewhere [19]. Note however that as the D13 has isospin 1/2 it does not contribute in the pn'^ channel and thus it does 
not affect the fit of the axial nucleon to Delta form factors carried out in Ref. [8]. 

For incoherent production on a nucleus we sum the nucleon cross section over all nucleons in the nucleus. For a 
charged current (CC) process, using the local density approximation, we aiTive for initial pion production (prior to any 
pion FSI) induced by a neutrino of momentum/energy |^| at 
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fC ) — y M^ + (^f ('■))^' being kp{r) = (3;r^pAr(r))'/^ and pAr(r) the local Fermi momentum and local density 

for nucleons of type A^. Besides 4>(|A;|) is the neutrino flux. d{vN — > l^N'n) is the cross section at the nucleon level 
modified by medium effects as discussed below. The above differential cross section is used in a simulation code to 
generate, at a given point r inside the nucleus and by neutrinos of a given energy, pions with a certain energy and 
momentum direction. 

The A properties are strongly modified in the nuclear medium [17, 20, 21, 22, 23, 24, 25, 26, 27] and since the direct 
A-pole contribution is the dominant one a more correct treatment is needed for production inside a nucleus. Following 
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FIGURE 2. Comparison of our model results (solid line) to ANL [5] and BNL [6] experimental data. Theoretical 68% confidence 
level bands are also shown. Data include a systematic error (20% for ANL and 10% for BNL data) that has been added in quadratures 
to the statistical published errors. Our theoretical results and ANL data include a W < 1.4 GeV cut in the final kN invariant mass. 



Ref. [23], we modify the A propagator in the A-pole term as 

1 1 



pl-Ml + iM^r^ V^ + Ma v?-MA + i(rP^"V2-imi:A)' 

with s ~ p\, r^™'' the free A width corrected by PauU blocking of the final nucleon, for which we take the expression 
in Eq.(15) of Ref. [27], and IoiEa the imaginary part of the A self-energy in the medium. The evaluation of Ea is done 
in Ref. [21] where the imaginary part is parameterized as 
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with po = 0.17fm^^. The terms in Ca2 and Ca^ are related to the two-body absorption WNN — > A^A^ and three-body 
absorption WNNN -^ A^A^A^ channels respectively. On the other hand the Cq term gives rise to a new WA^ — > Nn 
contribution inside the nuclear medium and thus it has to be taken into account beyond its role in modifying the A 
propagator This new contribution has to be added incoherently and we implement it in a approximate way by taking 

as amplitude square for this process the amplitude square of the direct A-pole contribution multiplied by y 12 • 

When coherent production on ^^C is possible we evaluate its contribution using our model in Ref. [17] but with the 
nucleon-to-Delta form factors as extracted in Ref [8]. 

As akeady mentioned, to evaluate FSI effects we follow Ref. [15] and we take into account P- and 5-wave pion 
absorption, and P-wave quasielastic scattering on a single nucleon. The P- wave interaction is mediated by the A 
resonance excitation. The different contributions to the imaginary part of its self-energy account for pion two- and 
three-nucleon absorption and quasielastic processes. The probabilities for the different processes are evaluated in 
nuclear matter as a function of the density and then the local density approximation prescription is used for its use in 
finite nuclei. After a quasielastic event, pions change momentum and may change its electric charge. The probability 
for charge exchange and the final momentum distribution after a quasielastic interaction are given in Ref. [15]. That 
information is used in the simulation program to generate the pion resulting from such a collision. Besides, in between 
collisions we assume the pions propagate in straight lines. All the details can be found in Ref. [15]. 



RESULTS AND COMPARISON WITH MINIBOONE DATA 



Here we show part of the results we have obtained. A more complete discussion of the relevance of the different 
contributions will be given in Ref. [19]. In the left panel of Fig. 3 we compare our results for n^ production in a CC 
process with MiniBooNE data. We take into account the contribution on ^^C and the two hydrogens. There is also a 
small coherent contribution on '^C Our total result is below data for neutrino energies above 0.9 GeV. The agreement 
improves if we do not take into account FSI of the pion. A similar result (see right panel of Fig. 3) is obtained for a 
final n'^. 
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FIGURE 3. l;r total production cross section for v^ CC interaction in mineral oil. Left Panel: results for a final n^ . Right panel: 
Results for a final ;r . Solid line: Total contribution. Double-dotted dashed line: Model prediction without FSI of the outgoing pion. 
Experimental data taken from Ref. [1]. 
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FIGURE 4. Differential ^ cross section for cliarged current In^ production by v^ in mineral oil. Captions as in Fig. 3. Data 
fromRef. [1]. 
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FIGURE 5. Differential ^^ (left panel) and ^^'^^^ (right panel) cross section for CC l^T*^ production by v^ in mineral oil. 
Captions as in Fig. 3. Data from Ref. [2]. 



In Fig. 4 we compare the differential J^ cross section for CC l;r+ production by v^. We have taken into account 
the neutrino flux in Ref. [1] to produce our results. We underestimate data for Tj^ above 0.15 GeV. This is an effect of 
FSI of the final pion which above those kinetic energies accounts for a sizable pion absorption driven by the A(1232). 
Neglecting FSI we get a good agreement with data. 

In Fig. 5 we show v^ differential jj- and j^'^q cross sections for CC Itt" production. For that we use the neutrino 

flux reported in Ref. [2] that extends from 2 GeV down to 0.5 GeV neutrino energy. Our results for jj- evaluated 
without FSI on the final pion agree better with data for pion momentum above 0.2 GeV/c. As a result of FSI, the 
agreement improves below 0.2GeV/c, but our model produces too few pions in the momentum region from 0.22 to 
0.55 GeV/c. The angular distribution shows those missing pions mainly go in the forward direction. 

In Fig. 6 we show results for NC production induced by neutrinos that we compare with data by the MiniBooNE 
Collaboration in Ref. [3]. We use the Vu flux reported by MiniBooNE. Our results for 4^ without FSI agree nicely 



with data, while our full model results show a depletion in the 0.25 
with data is nevertheless better than in the CC case. The differential 
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' 0.5GeV/c momentum region. The agreement 
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cross section is shown in the right panel 
of Fig. 6. Once again our results without FSI interaction of the final pion show a good agreement with experimental 
measurements. As for our full results, a clear deficit is seen in the forward direction but the agreement, as it was the 
case for the ^- differential cross section, is better than in the corresponding CC reaction. We obtain similar results 
for NC production induced by antineutrinos, see Ref. [19]. 

Our results both for CC and NC processes are in good agreement with the calculations in Refs. [10, 12]. As it is the 
case there, we also find a better agreement with data if FSI is ignored. The introduction of a formation time/zone, as 
done in Ref. [9], for pion production and its later interactions in the medium will decrease the effect of FSI and the 
agreement with data will improve. On the other hand, in Ref. [12] it is shown that the same FSI model applied to pion 
photoproduction on a nucleus is able to give a fair reproduction of experimental data. In Ref. [19] we also show that 
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FIGURE 6. Differential j^ (left panel) and j^^q (right panel) cross sections per nucleon for NC 1 K production by v^ in 
mineral oil. Captions as in Fig. 3. Data from Ref. [3]. 



our FSI model gives a fair reproduction of pion photoproduction in nuclei so that it is not clear to us what are the cause 
for disagreement in the neutrino induced reactions. 
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